BaVSa probed by V L edge X-ray absorption spectroscopy 
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Polarization dependent vanadium L edge X-ray absorption spectra of BaVSa single crystals are 
measured in the four phases of the compound. The difference between signals with the polarization 
E_Lc and E||c (linear dichroism) changes with temperature. Besides increasing intensity of one of the 
maxima, a new structure appears in the pre-edge region below the metal-insulator transition. More 
careful examination brings to light that the changes start already with pretransitional charge density 
wave fluctuations. Simple symmetry analysis suggests that the effect is related to rearrangements 
in Eg and Aig states, and is compatible with the formation of four inequivalent V sites along the 
V-S chain. 

PACS numbers: 78.40.-q, 71.30.+h, 71.27.+a, 71.45.Lr, 71.70.Ch, 71.45.-d 



I. INTRODUCTION 

Complex materials often display rich phase diagrams 
where different types of instabilities lead to unusual prop- 
erties related to various ordering phenomena. Under- 
standing the interplay between lattice, charge, spin and 
orbital degrees of freedom is one of most challenging as- 
pects of correlated electron systems. 

Barium vanadium sulfide, BaVSa, is a particularly in- 
triguing material because its structural, electrical and 
magnetic properties point to a combination of quasi-lD 
and 3D behavior. Its room temperature structure con- 
sists of hexagonal packing of quasi ID-chains of face- 
sharing VSg octahedra directed along the c axis 
The ratio of inter-chain to intra-chain V-V distances be- 
ing greater than 2, a high conductivity anisotropy ra- 
tio ((Tc/(Ta) can be expected. Surprisingly, ado a is only 
about three 0. Lowering the temperature induces a 
zig-zag deformation of the V chains at Tg = 240 K 
and reduces the crystal symmetry from hexagonal to or- 
thorhombic fl'] , but does not significantly change its con- 
ductive and magnetic properties. In both phases, it is 
a paramagnetic metal. Below Tmi = 69 K it becomes 
a paramagnetic insulator (PI). The metal-insulator (MI) 
transition is driven by a Peierls instability jsl, accompa- 
nied by a change to a monoclinic structure . At T/v = 
31 K, it undergoes a third phase transition to the lowest 
temperature (ground) state which has an incommensu- 
rate antiferromagnetic ( AF) order in the (a,b) plane . 
Recent resonant magnetic x-ray scattering reveal addi- 
tional incommensurability along the c axis in the AFI 
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ground state, while the supporting time-dependent den- 
sity functional theory simulations indicate that the spins 
lie within the (a,b) plane and are polarized along the 
monoclinic a axis [6|. 

Diffuse X-ray scattering experiments [3] have demon- 
strated that the MI transition is driven by a ID insta- 
bility of the electron gas. It leads to the formation of 
charge density waves (CDW) at the critical wave vector 
= 2fcF, preceded by pretransitional fluctuations. But 
the physical properties of BaVSs are not those of a stan- 
dard quasi-lD system. First, the wave vector q,.— 0.5 c* 
(c* stands for the reciprocal lattice vector related to c) 
indicates that only one of the two d electrons per unit 
cell (uc) participates directly in the CDW. On the other 
hand, instead of a typical Pauli behavior of a ID electron 
gas in the metallic phase, it has a Curie- Weiss behavior 
with an effective magnetic moment of /i = 1.2 /i^ cor- 
responding to approximately one localized spin per two 
vanadium sites All these properties point to an un- 
usual situation with one electron per uc delocalized along 
the c* direction driving the CDW formation. The other 
is localized and responsible for the magnetic behavior. 

Local density approximation (LDA) calculations indi- 
cate an interplay between two different types of electron 
states at the Fermi level: two narrow Eg bands, and 
one dispersive band with mainly A\g character extend- 
ing along the c* direction i. But the LDA filling of 
the A\g band is almost complete (n^i ~ 1-90 per uc), 
which overestimates the itinerant character, and is in- 
compatible with the above-mentioned experimental find- 
ings. Dynamical mean-field theory (DMFT) calculations 
show that strong enough charge correlation and exchange 
effects can bring ^i^i, down to one per uc [9]. 

Both types of electron states, delocalized (charge de- 
grees of freedom), and localized (magnetic properties), 
are modified at the MI transition. A charge gap of about 
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hPM 




Tjv= 31 K 


Tmi= 69 K Ts= 240 K 




structure 


monoclinic 


monoclinic 
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hexagonal 


space group 
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Cmc2i 




point group 
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m (Cih) 


mm2 {C2v) 


3m (Dad) 




4(4) 


4(4) 


2(1) 


2(1) 


magnetism 


AF 


para 


para 


para 


transport 


insulator 


insulator 


metal 


metal 



TABLE 1: Physical properties of BaVSa in its four phases: 
Anti- Ferromagnetic Insulator (AFI), monoclinic Paramag- 
netic Insulator (mPI), orthorhombic (oPM), and hexagonal 
(hPM) Paramagnetic Metal. Tc stands for critical tempera- 
ture, and Ny (Ni„eg.) the number of (inequivalent) V atoms 
per unit cell. 

50 meV and a spin gap of 10-20 meV [l^ are opened 
at TMI^ even if the spin degrees of freedom quench only 
in the AFI phase. But the mechanism of their interplay 
is not yet completely understood. 

We have explored four BaVSa phases in single crys- 
tals by polarization dependent x-ray absorption spec- 
troscopy (XAS) at the V L edges. After describing the 
experimental means for obtaining very high resolution 
contamination-free data we present the results and then 
discuss how the polarization dependence of the XAS data 
provides new insight into the subtle electronic structure 
changes that take place across a temperature range that 
covers all the phases of BaVSa. 

II. EXPERIMENT 

Single crystals of BaVSs were grown by the tellurium 
flux method [ll|. Two samples from different batches 
were cleaved in-situ along a plane parallel to the c-axis. 
After the cleaving at a pressure of 10~^ mbar, they were 
transferred via a fast insertion lock to the analysis cham- 
ber at the base pressure of 4 x 10~^° mbar. The experi- 
ments were performed using the helical undulator Dragon 
beam line ID08 at the European Synchrotron Radiation 
Facility (ESRF). At the V L3 edge (« 514 eV) the res- 
olution was w 100 mcV. The incident light was normal 
to the sample surface and linearly polarized either par- 
allel to the sample c-axis, or perpendicular to it. We 
periodically checked for energy shifts coming from the 
monochromator. The stability of the photon energy is 
estimated to be better than 40 meV. The V L edge XAS 
signal was measured in the total electron yield mode us- 
ing the drain current from the sample holder. 

III. RESULTS 

V L2.3 XAS results for BaVSs at room temperature 
for the two polarizations of the incoming light, E_Lc and 
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FIG. 1: (Color online) Paramagnetic Metallic phase (T — 
300 K) XAS spectra of the in-situ cleaved sample for the two 
polarizations of the incoming light: E_Lc (red dots) and E||c 
(black dots). The linear dichroism (LD = Ix-I||) is also shown 
in the same graph (green line). The inset shows the spectra 
of a sample cleaved in air where the oxygen contamination is 
clearly visible above 530 eV. Extra features are observed at 
the V L3 edge and the L2 peak is broadened. 



E||c are presented in Fig. [T] The spectra are normalized 
so as to obtain the same difference in amplitude below 
the L3 edge and at 534 eV. They are separated by the 2p 
spin orbit coupling, into L3 {hiy = 512-519 eV) and L2 
{hi' = 519-528 eV) regions. The V L3 and L2 maxima are 
at 516.2 eV and 523.2 eV, respectively, i.e. 7 eV apart. 
Their overall shape is closer to the spectra of vanadium 
metal [l^ than to the spectra of vanadium oxides that 
tend to be more structured (see for example VO2 [l^ or 

V2O3 [M])- 

The O K XAS signal lies in the 530-545 eV region. By 
comparing with the inset in Fig. [T] it can be seen that the 
surface contamination is negligibly small. In particular 
the hump at 540 eV is very weak while the feature at 
531 eV may even be intrinsic to the BaVS3 spectrum. A 
similar feature is present in elemental vanadium, where 
it has been attributed to a van Hove singularity at the 
boundary of the V metal Brillouin zone [15l] . The inset 
shows clearly that the oxygen contamination affects the 
V L edge. Additional structures show up at the high 
energy side of the V L3 peak, and the V L2 structure 
is broadened. In this sense, at least the surface of the 
slightly Ti-doped sample previously studied [ll] was con- 
taminated by oxygen. 

The V L2,3 maxima for the E||c polarization are at 
slightly higher photon energy than for the E_Lc polariza- 
tion. The shift is about 80 meV in the hPM phase. The 
peak position for V L3 with E±c excitation is 516.15 eV. 

In order to show more clearly the polarization depen- 
dence of the XAS intensity, we subtracted /|| (E||c spec- 
tra) from I± (E_Lc spectra). This linear dichroism (LD) 
curve is presented in Fig. [T] for the hPM phase. The 
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FIG. 2: (Color online) LD as a function of temperature. 
Red/light-green/dark-green/blue/black curves correspond re- 
spectively to the hPM/oPM-HT/oPM-LT/mPI/AFI temper- 
ature ranges. 
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FIG. 3: (Color online) LD changes with temperature in the 
pre-edge region. The spectra are shifted vertically for clarity. 
Inset shows details of the temperature changes in the oPM-LT 
range. 



change from positive to negative LD simply reflects the 
shift in the main XAS peak according to the polarization 
of the incident photons. This derivative-like shape hardly 
changes except for differences in amplitude as a function 
of the temperature or phase (Fig. [2]). 

When cooling to below Ts = 240 K (hPM to oPM), 
the shift increases to 110 meV. As a simplification we 
show the average spectrum taken at 300 K, 275 K and 
250 K to represent the hPM phase, what we term the high 
temperature orthorhombic phase (oPM-HT) is taken as 
the average of spectra recorded at 225 K, 200 K, 175 K, 
and 150 K, and the low temperature orthorhombic phase 
(oPM-LT) spectrum is the average of data taken at 
125 K, 100 K and 75 K. The mPI spectrum was recorded 
at 50 K and the AFI spectrum at 25 K. The shape of 
the main V L3 XAS peak remains practically unchanged 
but the amplitude of the LD signal at 515.8 eV changes 
reflecting very small variations in the relative intensity of 
the /|| and I± maxima. Another peak in the LD is ob- 
served at 514.2 eV, coinciding with the shoulder on the 
XAS data. Its amplitude is positive and remains practi- 
cally unchanged across the range of measurements. 

An essentially equivalent situation is observed for the 
L2 edge. Here structures are not well resolved. This can 
be explained by the difference in lifetime broadening of 
the two edges, which is 0.78 eV (0.28 eV) for the V L2 
(L3) [17]. 

A detailed examination of the low energy part of the 
spectra, brings to light structure at 512.8 eV. The LD 
data are shown in Fig. [31 As the temperature is lowered 
to the oPM-LT temperature range the intensity is in- 
creased at 513.2 eV, where a relatively strong extra peak 
is clearly present in both the mPI and AFI phases. 

It should be underlined that we have observed the same 
low temperature behavior for two samples from different 
preparations. Fig. 2] shows the LD for the sample from 
the second batch demonstrating that the data are highly 
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FIG. 4: (Color online) LD changes with temperature in the 
pre-edge region for the sample from another preparation is 
essentially the same. 



reproducible. 



IV. DISCUSSION 

An elaborate theoretical interpretation of the BaVSa 
XAS spectrum is not presently available because of its 
structural complexity. We will therefore discuss our data 
in terms of a simplified analysis of orbital symmetries 
to help assess what information can be extracted from 
the XAS polarization dependence, along with our present 
knowledge of the physical properties of the material and 
the insight that can be gained from recent LDA-DMFT 
calculations fisj . 

First we need to calculate the polarization dependence 
of the 2p to 3d dipole transition in the material to assess 
the polarization effects. 
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Due to the inclination of the V-S octahedra in BaVSs, 
it is convenient to present its V 3d states as a linear 
combination of cubic d orbitals (Eq. [T] and Fig. [5|): 



Aig = dz 



En 



En 



dx2_y2 



2 ^x^—y^ 



(1) 



The 2p — > 3d transition probabilities, according to dipolar 
selection rules, are given by 



P, = ^|(KlE-r|, 



(2) 



where pi represents three V 2p orbitals and d can be any 
of the five V 3d orbitals: Aig, Egi, Eg2, egi, eg2- E • r is 
the dipolar transition operator. 




FIG. 5: (Color online) Five BaVSs V 3d orbitals correspond- 
ing to Eq. [1] two Sg (upper panels) and three t2g (lower pan- 
els). They are inserted in the orthorhombic phase deformed 
octahedra, with two inequivalent sulfur atoms. Si (orange) 
and S2 (yellow). Axes (x,y,z) correspond to (a,b,c) crystallo- 
graphic axes in the orthorhombic phase. 



The polarization dependent transition probabilities, 
for E±c and for E||c, are the following: 

- Transitions to the higher energy Cg states are not po- 
larization dependent. For both polarizations their prob- 
ability is P^^ = pI = 0.4. 

- Transitions to the Aig states are strongly polarization 
dependent, Pj[^^ = 0.1 and p|^_^ = 0.4. 

- Transitions to Eg states are also polarization dependent, 
P^ = 0.5 and pj, = 0.2. Thus the total polarization 
dependent absorption intensity as a function of energy 
can be written as: 
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FIG. 6: Local spectral functions from LDA+DMFT calcu- 
lations of the t2g states in the oPM phase where all V sites 
are equivalent (upper panel) and mean value of local spectral 
functions for four inequivalent V-sites for the t2g states in the 
mPI phase (lower panel). The LDA density of Cg states shown 
in both panels is the one calculated for the oPM phase, as the 
calculation for the mPI phase is not available. Fermi level 
energy is set to 0. LDA and DMFT data have been adapted 
from the Figs. 2, 7 and 12 of RefjU. 



I±,\\ie) cx nA,,ie)P^;l + K,, (e) + n^^, (e)]P4'" 

-fK,(e)-|-ne,,(e)]P,-^-ll, (3) 

where nd(e) are unoccupied spectral function weights 
at a particular final (core-excited) state energy e. 

As the transitions to the eg states are equiprobablc for 
both polarizations, the linear dichroism should be zero for 
excitations to states with this symmetry. The difference 
in the transition probabilities to Aig is equal to that of 
each individual Eg orbital but of opposite sign. Thus 
an increase in the population of Aig states will tend to 
reduce the LD or give a negative signal while an increase 
in the Eg population will tend to increase the LD when 
defined as (/_l — /||)(e). 

Summarizing our observations, little difference in the 
polarization dependence can be detected across the 
Ts = 240 K hPM to oPM phase transition despite the 
symmetry break which lifts the degeneracy of the two Eg 
and two Cg states as shown in Fig. [6] (upper panel). Low- 
ering the temperature below 150 K (oPM-LT) induces 
an increase in the LD amplitude at 515.8 eV, while the 
intensity starts to build up at 513.2 eV, where a distinct 
peak appears below Tmi — 69 K. When the temperature 
is further reduced below the mPI to AFI {T^ = 30 K) 
phase transition the intensity at 515.8 eV again drops 
and the 513.2 eV peak intensity strengthens. This could 
mean that spin ordering has the effect of transfering one 
part of unoccipied Eg states to lower energy. 

The increase in LD amplitude at 513.2 eV and 515.8 eV 
across the oPM-HT to oPM-LT temperature range might 
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be interpreted as an increase of the Eg weight in the cor- 
responding small energy range. As the magnetic prop- 
erties are related to the Eg states, this interpretation is 
in line with the energy dependence of the ground state 
magnetic reflection at the (0.226 0.226 0.033) wave vector 
which shows two V L3 peaks at similar positions 6]. On 
the other hand, in the oPM-LT temperature range, there 
is a formation of CDWs in chains, as shown by pretransi- 
tional diffuse lines in x-ray diagrams Q . It is a dynamical 
effect, with a typical lifetime of 10""'^^ s [l^, which can 
induce local structure modifications via electron-phonon 
interaction. As the time scale of XAS is determined by 
the V L3 edge core-hole-lifetime of 10~^^ s, it can probe 
short-time modifications of unoccupied states related to 
CDW fluctuations. Similar dynamical local structural 
change is observed at the Ni K-edge NEXAFS of PrNiOa 
[20j system. 

Over the oPM-HT temperature range all V sites are 
equivalent by symmetry, while in the mPI phase there 
are four inequivalent V-sites in the monoclinic unit cell 
[21I I . Resonant diffraction measurements at the V K edge 
pointed to the possibility of an orbital ordering of AOEO 
type (Aig, equal, Egi, equal) accompanying the CDW 
formation [2^]. Collective intersite orbital excitations 
measured in the dielectric signal support this hypothesis 
[2^. DMFT calculations confirm that the four inequiva- 
lent V-sites have different spectral functions but propose 
rather an EEOO occupation pattern [l^. In terms of 
unoccupied states, two sites have their Eg2 local spec- 
tral function maximum closer to E^? and the other pair 
have the maximum pushed to higher energy, giving rise 
to a broad double peak when the mean value for four 
inequivalent V-sites is presented as in the Fig. |6] (lower 
panel). 

From these considerations we suppose that, below 
150 K, the increased LD intensity at 513.2 eV and 
515.8 eV is related to the strong high energy transfer 
of a part of unoccupied Eg states. This shift is itself re- 
lated to the formation of four inequivalent V-sites which 
accompanyies CDW formation in the pretransitional fiuc- 



tuation regime. Possible core-hole effects are difficult to 
estimate without proper theoretical modeling. It is how- 
ever reasonable to suppose that difference in core hole 
interaction across the phases is negligibly small. 



V. CONCLUSION 

BaVSa single crystal V L edge x-ray linear dichroism 
(LD) spectra show changes with temperature. Below 
150 K there is a clear LD intensity increase at 513.2 eV 
and 515.8 eV. Simple symmetry analysis suggests the ef- 
fect is related to rearrangements in Eg and Aig states. 
Two V L3 peaks at similar energies are observed in the 
energy scan of the ground state (0.226 0.226 0.033) re- 
flection measured by the resonant magnetic x-ray scat- 
tering . This confirms that the observed increase of the 
LD intensity is related to magnetically active Eg states. 
DMFT calculations [18j predict rearrangements of Eg 
states in the mPI phase, together with CDW-related for- 
mation of four inequivalent vanadium sites with different 
Eg and Aig spectral functions. Our data show that the 
Eg rearrangements start at higher temperatures with the 
formation of CDWs in chains, affecting Aig states. We 
expect these experimental results to help ongoing efforts 
to model this complex material as they highlight the need 
to consider Aig and Eg states as a whole to understand 
the physical properties of this unusual compound. 
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